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C/^ I Abstract. The electrovacuum around a rotating massive body with electric charge density is 

. described by its multipole moments (mass moments, mass-current moments, electric moments, 

fS| ' and magnetic moments). A small uncharged test particle orbiting around such a body moves on 

\ geodesies if gravitational radiation is ignored. The waves emitted by the small body carry informa- 

T-H . tion about the geometry of the central object, and hence, in principle, we can infer all its multipole 

' moments. Due to its axisymmetry the source is characterized now by four families of scalar mul- 

, tipole moments: its mass moments M;, its mass-current moments Si, its electrical moments £/ and 

■ its magnetic moments Hi, where / = 0, 1,2, . . .. Four measurable quantities, the energy emitted by 

gravitational waves per logarithmic interval of frequency, the precession of the periastron (assuming 

almost circular orbits), the precession of the orbital plane (assuming almost equatorial orbits), and 

' the number of cycles emitted per logarithmic interval of frequency, are presented as power series of 

\ the newtonian orbital velocity of the test body. The power series coefficients are simple polynomials 

— ■ of the various moments. 
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: INTRODUCTION 

\^ • 

c3 ■ The geometry around a massive object is uniquely characterized byo the muhipole mo- 

ments [L2».3„ 4] of the central object, as in newtonian gravity. Ryan |5] showed how 
we could reveal the spacetime geometry of a central axisymmetric body, by measuring 
a few physical quantities that are related to the kinematics of a hypothetical test body 
that is orbiting around the central object while it is emitting gravitational waves. This 
geometry-mapping could be achieved in principle by an advanced network of gravita- 
tional wave detectors with high sensitivity. Of course, to determine all moments we need 
infinitely sensitive detectors. However, knowing a few moments of the central body, we 
could put restrictions on the various models that are assumed to describe the interior of 
the central body. 

Since the electromagnetic moments of the central object have their own contribution 
on the geometry around the central object, we could generalize the work of Ryan to yield 
the relation of all moments to the physical quantities used by Ryan. 

In our paper we have considered binaries consisting of (i) a central object is assumed 
to be stationary and axisymmetric, and characterized by reflection symmetry with re- 



spect to its equatorial plane, (ii) a test particle with no charge that is orbiting around 
the first one on a nearly circular and equatorial geodesic orbit. The axisymmetry of the 
central object ensures that the spacetime geometry could be described by scalar mo- 
ments. Since we are dealing with electrovaccuum solutions 4 sets of scalar moments are 
needed: Its mass moments Mo,M2,M4, . . ., its mass-current moments Si, S3, . . ., its elec- 
tric charge moments Eq,E\,E2, . ■ and its magnetic moments Hq, Hi,H2,..., where e.g., 
the Mq = M moment is the mass of the object, 5i is its angular momentum, Eo = Eis its 
charge, and Hi is its magnetic dipole moment. In every set of moments, each moment 
appears in steps of 2, and this holds good for the electromagnetic moments as well, due 
to the reflection symmetry of the metric itself (cf., [2]). The fact that all electromag- 
netic moments are written above, instead of every second one, is related to the fact that 
two kind of objects could have axisymmetric metrics, either the ones with even electric 
moments and odd magnetic moments, or the other way around. 

The two precessing frequencies of the orbits (the periastron precession and the orbital 
precession) are written as functions of the moments of the metric on which the particle 
is moving. Moreover, although the orbit of the test body is assumed to be a geodesic one, 
we are going to use the quadrupole formula to infer the energy emitted by the binary. 
The evolution of the orbit will then be used to obtain alternative information about the 
moments of the geometry. 

Throughout the paper units are chosen so that G = c = 1 . 



THE KINEMATICS OF THE TEST PARTICLE 

If we perturb a circular equatorial orbit on a fixed axisymmetric background then two 
new frequencies arise, the one connected to the periastron precession and the one 
connected to the orbital precession, that are related to the circular orbit frequency and 
the metric itself by |H] 
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where a stands for p (periastron precession), or z (orbital precession). Actually, the 
frequencies written above correspond to the difference between the orbital frequency 
and the frequency of perturbations in p, or z, since these differences will show up in 
gravitational radiation emitted by the binary as a modulating frequency. 



The energy per unit test-body mass for an equatorial circular orbit in an axially 
symmetric spacetime is 



u ~ I ' 



and the specific energy released as gravitational radiation per logarithmic interval of 
frequency is 

^ = (3) 

The expression above assumes that all the energy lost from the test body has been 
emitted at infinity as radiation. 

The number of gravitational-wave cycles spent in a logarithmic interval of frequency, 
the quantity that could be most accurately measured by a gravitational wave detector, is 

"^wave/ 

where dE^f^aMe/dt is the gravitational-wave luminosity, which is assumed to be exactly 
the rate of energy loss of the orbiting test body. The main contribution of dE^a^y^/dt 
comes from the mass quadrupole radiative moment of the binary. Therefore, when we 
give the connection between the various orders of v = {MQ)^/^ (the orbital velocity) in 
the power series expansion of /SN, and the moments of spacetime, only the highest order 
moments enter our expressions, since lower order moments contribute through higher 
than quadrupole radiative moments as well in a rather complicated and unclear way. 
However, at least the first 4 power parameters of /SN could be written explicitly in terms 
of the multipole moments. 



THE METRIC AS A POWER SERIES OF MOMENTS 

To get expressions for these 3 measurable quantities, that are straightforwardly con- 
nected to the moments of the central object, one has to reexpress the metric functions in 
terms of all moments, as well as the radius p and the orbital frequency of the test body 

In our paper, we consider only stationary axisymmetric objects that are reflection 
symmetric with respect to their equatorial plane. The metric of such a central object 
alone could be written in p,z, ^) coordinates, in the form of Papapetrou metric 

ds^ = -F{dt - CO d<^f + J [e^^{dp^ + dz^) + p^d^^] , (5) 

where F, co, and 7 are the three functions that fully determine a specific metric. These 
are functions of p and \z\ only, due to axisymmetry and reflection symmetry. Einstein's 
equations in vacuum guarantee that once F and CO are given, 7 can be easily computed 
(see Once we incorporate an electromagnetic field in the vacuum around the com- 
pact object, the source of which is the compact object itself, which allows spacetime to 



have the same symmetries, the metric above still describes the electrovacuum spacetime, 
but now the metric and the electromagnetic field should satisfy the Einstein-Maxwell 
equations. In order to fully compute the metric functions one more complex function, 
4>, which is related to the electromagnetic field, is necessary. F , CO, and 4> themselves can 
be determined by solving the so-called Ernst equations ||E0]> which are the Einstein- 
Maxwell equations written in a different form. It is a system of non-linear complex 
differential equations of second order: 

(9t(f?)-F|4>|2)vV = (V^ + 24>*V4>)-Vff, (6) 
{3i{<ff) + \^\^)V^^ = (Vcf + 24)*V4>)-V4>, (7) 

where V denotes the gradient in a cartesian 3D space (p,z, ^) and . .),3(. . .), here 
and henceforth, denote the real and imaginary part, respectively, of the complex function 
in parentheses, while an asterisc * denotes complex conjugate. The third metric function, 
7, is then easily computed by integrating the partial derivatives dy/dp, dy/dz, which 
are given as functions of derivatives of all other functions (see ifioll '). 

From the two complex functions, <S and 4>, we could construct the metric functions as 

S = (F-|4>p)-fi(p, (8) 

where ^ is related with gt^ through 

(9) 
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Note that there is a sign difference in Eq. (22) of Jsl], which has been corrected in a later 
paper of Rva n lllh . and comes from an odd convention of CO used by Ernst (see relevant 
comment of WM ). 

Now, instead of <S' and 4>, one could use two new complex functions t, and q, that 
play the role of gravitational potential and Coulomb potential respectively, and are more 
directly connected to the mass and electromagnetic moments of the central body. These 
potentials are related to the Ernst functions by 



s = ^/LZl \ (10) 
^ = ^^^^ 



and can be written as power series expansions at infinity 
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and cijj.bij are coefficients tliat vanish when z is odd. This reflects the analyticity of the 
potentials on the z-axis. The tilded quantities, here and henceforth, are the conformally 
transformed ones, which are essential for calculating the moments (see [1^)- 

Due to Ernst equations (I6l7l) the above power expansion coefficients, atj and bij, are 
interrelated through complicated recursive relations (see ifujl V Essentially, these rela- 
tions are simply an algebraic version of Einstein-Maxwell equations for the coefficients 
of the power expansion of the metric and the electromagnetic field tensor. The recursive 
relations could be used to build the whole power series of ^ and q from a full knowledge 
of the metric on the axis of symmetry 

oo oo 

I (P = 0) = £ m-z\ q{p = 0) = £ . (14) 
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In 111311 a method of calculating the complex multipole moments of the central object 
in terms of the m/s and ^/s is presented. In brief, the gravitational moments are given 
by 

^ St\ (15) 
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where Sa are given by recursive formulas analyzed in 

The mass moments M„ and the mass-current moments 5„ are related to by 

Pn=Mn + iSn, (16) 

whereas the electric moments En and the magnetic moments H„ are related to Q„ by 

Qn = En + iHn. (17) 

Since this algorithm can be used to evaluate the moments in terms of the m, 's and qi's, 
one can invert these relations and express the m,'s and qis in terms of the moments: 

mn = aon = Mn + z5„ + LOM, 

qn = bon = En + iHn + LOM, (18) 

where "LOM" stands for lower order multipole moments of any type. Thus, we can use 
the recursive relations mentioned above to evaluate the ajj and bij coefficients in terms 
of the moments. 

Finally, following the procedure presented in the beginning of this section we can 
express the metric functions, and their first and second derivatives as power series of p 
and z with coefficients that are simple algebraic functions of the moments of the massive 
body. Since in our study we have confined the motion of the test particle on the equatorial 
plane, we actually need to compute everything at z = which makes calculations far 
simpler than what they seem. 

One could argument (see Sec. IIC of [10]) that the reflection symmetry assumed for 
the metric is consistent with a set of even electric and odd magnetic moments or odd 
electric and even magnetic moments. 



EXPRESSIONS RELATING THE MEASURABLE QUANTITIES 

WITH MOMENTS 



Combining the formulae that are given in Sec. 2 we can express the three measurable 
quantities as power series of v = (MQ.) with coefficients that have explicit dependence 
on all four types of moments. The choice of v as a dimensionless parameter to expand 
all physical quantities is warranted from the fact that the inspiral phase of a binary, 
the best exploitable part in gravitational- wave analysis ifTSl . involves comparatively low 
magnitudes of v. All measurable quantities have been transformed to a dimensionless 
form as well, for example by dividing the two frequencies Q.p, by the orbital 
frequency Q.. 

We need also a power series expansion of p with respect to v, or equivalently Q.. Thus 
we have to invert the function 0.{p), at least as a power expansion. From an elementary 
analysis of circular geodesies on the equatorial plane (see ||5|]) we know that 



a = ^ . (19) 

g<i>(i),p 

In the following part of this section we explain the algorithm that one should follow, in 
order to obtain the power series for Q.p/Q., Q.^/ Q., and AN. One starts with a power series 

of I and q of the form given by Eq. dT^ . Since no higher than second derivatives of the 
metric functions with respect to z are necessary, one should keep at/s and bi/s with 
< J < 2, and as many values of / as one needs to carry the power series expansion of 
the measurable quantities at a desirable order. All quantities that are expressed as power 
series of p and z, are evaluated at f = z = at the end, and thus, all expressions are finally 
power series of p = 1/p, due to Eq. (fT3l) . Although the ciij and bij are polynomials of 
various moments, from the practical point of view it is preferable to keep them as they 
are, and replace them by their moments dependence only at the final expressions. Then 
from ^ and q we construct S', 4>, and F, cp (cf., Eqs. (I8I10I) '). These are sufficient to build 
all metric functions. Next, following the procedure described above, we expand Q. as a 
power series of 1/p, by virtue of Eq. dT^ . This series is inverted and in this way we 
obtain 1 /p as a power series of Q., which then can easily be turned into a power series 
of the dimensionless parameter v. 

Now, the power series representing 1/p will replace all 1/p terms appearing at the 
expansions of the metric, its derivatives, and all other physical quantities depending on 
them. Finally, one has to rewrite the atj and bij terms appearing at the coefficients of all 
these power series as polynomials of the various moments. The recursive relations that 
relate all ajj and bjj with mk = aok and q^ = bok, which are directly related to the scalar 
moments of spacetime through Eqs. (24,25) of Ref. [13]. 

The algorithm described in the previous two paragraphs has been carried out with 
Mathematica, and has been checked for the following two subcases: (i) When all elec- 
tromagnetic fields are turned off, by erasing all electromagnetic moments (Ei = Hi = 0), 
our expressions for Qp, Q.y, AN are identical to the ones computed by Ryan [5]. (ii) 
For the Kerr-Newman metric it is quite easy to compute Q., Q.p, Q.^, and AE / fi for 



a quasi-equatorial, quasi-circular orbit. The expressions we obtain are identical to the 
ones obtained directly from the Kerr-Newman metric. 
The power series expansion for Q.p, Q.^, take the following form: 

^ = Ii?nV«, (20) 

§ = t^nV", (21) 

n=3 

while the corresponding coefficients, up to 5th order for the two frequencies, for the two 
distinct electromagnetic cases (es) and (ms) are 
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The fact that we have two new sets of moments (the electromagnetic ones) with 
respect to Ryan allows many more combinations of moments in high order terms. 
Actually, from a practical point of view these expansions are far more advanced than 
what will be used in gravitational wave data analysis in the near future. On the other 
hand the expressions above present an important feature: in every new order term a new 
moment shows up. This suggests that a very accurate observational estimation of the 
series could in principle reveal any moment. 

A glance at the corresponding terms of the two electromagnetic cases shows that each 
combination of moments for the (es) case is numerically equal to the corresponding com- 
bination for the (ms) case, if the electric and magnetic moments are interchanged. The 
sign though is the same for combinations of pure electric, or pure magnetic moments, 
but opposite for combinations of electric and magnetic moments. 

Finally, in order to express AN also as a power series of v, we need to expand 
dE\i/ave/dt as power series of v. As was explained in Sec. 2 we cannot work out the 



perturbative analysis of gravitational wave emission at a generic spacetime background; 
we can only obtain accurate expressions for JE'wave/'^^ up to after the leading or- 
der. Therefore, in the following formulae for p and dE^asQ/dt, we write simply the 
power series coefficients explicitly up to the fourth order. More specifically, in order 
to compute the power expansion of J^wave/^^ we add up all power series contribu- 
tions (quadrupole radiative moment, current quadrupole radiative moment, and Post- 
Newtonian corrections) which are the only ones that show up up to the 4th order. Thus, 
we yield 

^-Miiy['^ '''' 

where the coefficients for the two electromagnetic cases are given by the following 
polynomials of the moments 
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As in the rest three measurable quantities, the power expansion of AN is such that in 
every order term a new moment, which was not present in any lower order term, occurs. 
This is an indication that all moments can in principle be unambiguously extracted from 
accurate measurements of AN. 



COMMENTS 

As is shown in flsl the first generation of LIGO is not expected to be able to extract 
the first two moments (Si and M2) with high accuracy (~ 0.05 for the former and ~ 0.5 
for the latter one), by analyzing the phase of the waves. If we allow for electromagnetic 
fields as well, the corresponding monopole (which classically is expected to be very 



close to zero) will be measured with even higher accuracy than the other two mass 
moments, since the charge of the source (or the magnetic monopole in case of some 
exotic body) is present at even lower order, namely in the term, while the electric 
dipole, or the magnetic dipole, that first show up at the term will be measured with 
rather disappointing accuracy. On the other hand, analyzing the data of LISA leads to 
accuracies almost two orders of magnitude higher than the corresponding for LIGO. 
Thus, it seems quite promising that LISA will give us the opportunity to measure the first 
few moments, including the electromagnetic dipole moments, quite accurately. Also, 
the fact that in every new order in the power series of AA'^ a new moment appears is 
significant, since this means that in principle a unique set of moments arises from an 
accurate estimation of all power series terms. Actually, there are two possible sets of 
moments; one for each electromagnetic case, since we cannot a priori exclude one of 
them. We can only exclude one of the two sets on physical grounds, if only one of 
them leads to a physically reasonable classical object (for example a highly magnetized 
compact object is physically preferable to a compact object with a huge electric dipole). 
If we manage to measure a few lower moments, we can check if they are interrelated 
as in a Kerr-Newman metric llisl] . A positive outcome of such a test will be of support 
to the black-hole no-hair theorem in the case that the central object is a black-hole. The 
case of observational violation of the black-hole no-hair theorem could either mean that 
the central object is not a black-hole, or that the theorem does not hold. Of course to 
assume the latter an extra verification that the central object is indeed a black hole is 
necessary. 

While the phase of a gravitational wave is the quantity that can be most accurately 
measured, since a large number of cycles (a few thousand for LIGO and a few hundred 
thousand for LISA in case of binaries with high-ratio of masses) is sweeping up the sen- 
sitive part of the detectors, the two precession frequencies Q.p and Q.^, can in principle 
be measured if the detectors become more sensitive and templates that describe modu- 
lating waves are used lH^ . If this ever become possible one could use any of them to test 
the no-hair theorem. This would demand no more than the four lower order terms, since 
according to this theorem all moments depend on only three quantities (mass, angular 
momentum, and total charge). Actually, from measurements of modulating frequencies 
we could not at first determine which frequency corresponds to each precession. How- 
ever, the power expansions of the two frequencies begin at a different order, and thus we 
could discern them. Unfortunately, there are terms in these expansions that contain more 
than one first occurring moment (for example the first term of Q.p, 7?2, is a function of E 
ovH and M). However, expansions of and Q.y, if used simultaneously, along with the 
intrinsic dependence of v on M, could finally lead to the full determination of the mo- 
ments. Notice though that some corresponding terms in the series expansions of the two 
frequencies depend on the same set of multipole moments, like and Z3. Assuming 
that the observed system is sufficiently well described by our model for the binary, this 
multipole information could serve to determine these moments with better precision. 

If we manage to obtain information about the precessing frequencies or the evolution 
of the orbits of small bodies around compact stars by other observational means, we 
could measure the multipole moments (mass moments, mass-current moments, and 
electromagnetic moments) of the compact central object and thus put strict restrictions 
on the models that are used to describe the interior of these objects. 



Our analysis demonstrates that it will be possible to determine all types of multipole 
moments of the central object, from future gravitational wave measurements. Thus, apart 
from spacetime geometry, we could also determine the central body's electromagnetic 
fields. Although the data of LISA should be suitable for extracting such information 
with high accuracy, the assumptions of circular and equatorial orbit are not that realistic. 
From this point of view we consider our work as a step towards a more detailed analysis 
with not so restrictive assumptions. 
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